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Abstract. We develop a Hamilton- Jacobi theory for singular la- 
grangian systems in the Skinner-Rusk formalism. Comparisons 
with the Hamilton- Jacobi problem in the lagrangian and hamil- 
tonian settings are discussed. 
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1. Introduction 

The standard formulation of the Hamilton- Jacobi problem is to find 
a function S(t, q A ) (called the principal function) such that 

apt) = o, (i-i) 



Key words and phrases. Hamilton- Jacobi theory, presymplectic constraint 
algorithm. 
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where h = h(q ,Pa) is the hamiltonian function of the system. If we 
put S{t, q A ) = W{q A ) - tE, where E is a constant, then W satisfies 

h{ q \ d ^ K ) = E- (1.2) 

W is called the characteristic function. 

Equations (II. ip and (11. 2p are indistinctly referred as the Hamilton- 
Jacobi equation (see [Ij |2"| 125]). 

This theory works for classical mechanical systems, where the la- 
grangian function is usually the kinetic energy corresponding to a Rie- 
mannian metric on the configuration manifold minus a potential energy. 
This is the case of the so-called regular lagrangian systems, that have 
a well-defined hamiltonian counterpart. The theory has been recently 
reformulated in a geometrical setting (see [SHU [3]) that has permitted 
its extension to nonholomic mechanical systems [13] [15] , and even clas- 
sical field theories [TH [18] .The procedure is based on the comparison 
of the hamiltonian vector field Xh on the cotangent bundle T*Q and 
its projection onto Q via a closed 1-form 7 on Q; the result says that 
both vector fields are 7-related if and if the Hamilton- Jacobi equations 
d(h o 7) = holds. 

On the other hand, a Hamilton- Jacobi theory for singular lagrangian 
systems is far to be accomplished. There were several attempts ([20 ] [21] 
[22]). based on the homogeneization of the given lagrangian, which leads 
to a new lagrangian system with null energy such that it is possible to 
discuss the Hamilton- Jacobi equation for the constraints themselves. 
The main problem is that, due to the integrability condition for the 
resultant partial differential equation, one can only consider first class 
constraints. Therefore, the treatment of the cases when second class 
constraints appear should be developed by ad hoc arguments (as in 
[22] . for instance). Thus, in [20] and [21] the authors only discuss the 
case of primary constraints. 

A more modern discussion on this subject can be found in [3J fT2"] . 
but these authors only consider the case of primary constraints. More 
recently, in [T7] it is proposed a Hamilton- Jacobi theory for arbitrary 
singular systems that works even if the system exhibit secondary con- 
straints. The strategy is to apply the geometric procedure described 
above in combination with the constraint algoritm developed by M.J. 
Gotay and J.M. Nester [TJ [HI [10] and that geometrizes the well-known 
Dirac theory of constraints [5]. 

In the present paper we take a different approach, and consider the 
Skinner and Rusk setting to treat with singular lagrangians [2J1 [25] . 
Skinner and Rusk have considered a geometrized framework where the 
velocities and the momenta are independent coordinates. To do this, 
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they considered the dynamics on the Withney sum of TQ (the space 
of velocities) and T*Q (the phase space). 

Given a lagrangian function L : TQ — > R (singular or regular, no 
matter) one considers the bundle TQ@T*Q with canonical projections 
W\ ■ TQ © T*Q TQ and pr 2 : TQ © T*Q -> T*Q onto the first 
and second factors. We then define a function D : TQ © T*Q — > R 
by D(X p ,a p ) = a p (X p ) — L(X p ). In bundle coordinates (q A , v a ,pa), 
D is given by D(q A , v a ,Pa) — ^ A Pa — L(q A ,v A ), and it is sometimes 
refered as the Pontryagin hamiltonian or generalized energy (see |26j). 
We can also define a 2-form Q on TQ © T*Q by Q = pr^fig), where 
Qq denotes the canonical symplectic 2-form of T*Q. 

Then, one discuss the presymplectic system (TQ © T*Q, Q, dD) and 
obtain the corresponding sequence of constraint submanifolds, which, 
of course, have a close relation with those obtained by Gotay and Nester 
on the lagrangian and hamiltonian sides. It should be noticed that this 
algorithm includes the SODE condition just from the very beginning. 

We apply the Hamilton- Jacobi geometric procedure to this presym- 
plectic system and develop the corresponding Hamilton- Jacbi theory. 
The relation with the Hamilton- Jacobi problems on the lagrangian and 
hamiltonian sides are extensively discussed. 



2. Notation and background 

In this work all manifolds are assumed to be finite dimensional and 
C°° . Given a function /, the differential at a point p will be indistinctly 
denoted by d p f or df(p). 

We refer to [19] for a detailed description of lagrangian and hamil- 
tonian mechanical systems. 

Let Q be a differentiable manifold and denote by TQ and T*Q the 
tangent and cotangent bundles, and by tq : TQ — > Q and txq : T*Q — > 
Q the respective canonical projections on Q. 

We introduce two canonical structures on the tangent bundle of a 
manifold: the vertical endomorfism S, and the Liouville vector field A. 
In bundle coordinates, (q A ,v A ), they are respectively given by 

S = dq A ®^, 
A = v A 

Let now L : TQ — > R be a lagrangian on TQ; we can define the 
Poincare-Cartan 2-form and the energy function of L by 



Q L = -dO L , where 9 L = S*{dL), 
E L = A(L) - L, 
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which in local coordinates read as 

n L = dq A Ad-^ J , 

a 9L T . , 
E L = v q^a ~ L(q,v). 

We look for vector fields £ which simultaneously satisfy the equations 

i 6 n L = dE L (2.1) 

S^ = A. (2.2) 

If the lagrangian L is regular, that is, det( di ^g v B ) 7^ 0, then the form 
Ql is symplectic (Ql has maximal rank) and there exists a unique 
vector field £ on TQ which satisfies the equation (12. ip . This vector 
field automatically satisfies the SODE condition (12.21) . 

If the lagrangian is not regular, then Ql is no longer symplectic and 
equation (12.11) has no solution in general and even if there is a solution 
it is not necesary a SODE. Therefore for a singular lagrangian L, 
is a presymplectic form (that is, the rank is not maximal, althought, 
for simplicity, it is assumed that it is constant). 

We define the Legendre transformation associated to L as the map- 
ping 



FL : TQ — ► T*Q 

' dv A 



{q A ivA) ^ FL(q A ,v A ) = (q A ^(q A ,v A )). 



From a direct inspection in local coordinates we know that the Legendre 
transformation is a local diffeomorfism if and only if L is regular. 

We can apply the Gotay-Nester-Hinds algorithm of constraints, see 
El E], to the presymplectic system (TQ, Q L , dE L ) and hence we 
obtain a sequence of constraint submanifolds 



• • P k ^ ■ ■ ■ P 2 P 1 = TQ. 



Assume that the algorithm stabilizes at some step k, say Pk+i = Pk, 
which is called the final constraint submanifold, denoted by Pf = Pk- 

In this paper we will only consider almost regular lagrangians L : 
TQ ->■ R, that is: 

(i) Mi =Im(FL) is a submanifold of T*Q, and 

(ii) FL : TQ — > Im(FL) is a surjective submersion of connected 
fibers. 
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Under these assumptions, the energy E L is projected onto a function 
h x : Mi ->■ R such that hi o FL = E L 




Mi = lm(FL) 

Here FLi is the restriction of FL to its image, and ji : Mi — > T*Q is 
the canonical inclusion. 

Next, study the presymplectic system given by (M 1; Qi = jI^Iq, dhi), 
where fig is the canonical symplectic form on T*Q. Therefore, we con- 
sider the equation 

i Y Sli = dhi. (2.3) 

As above we can apply the presymplectic algorithm and we obtain a 
sequence of constraint submanifolds 



• •-.!/,, M, ■ y M : ^ > T*Q. 



It is obvious that 

FL(Pi) = M^ for any t, 

and, furthermore, the induced mappings 

FLi = FL\ P . :Pi^Mi 

are surjective submersions, for all %. 

Hence, both algorithms stabilizes at the same step, say k, and then 

FL(P f ) = M f , 

and 

FL f :P f ^M f 
is a surjective submersion (with the obvious notations). 
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The following diagram summarizes the above discussion. 



Pi = TQ 



FL 



T*Q 




Mi 



M 2 




M. 



where gi and jj denote the natural inclusions. 

The relation between equations (12. ip and (12. 3p is given by the fol- 
lowing theorem. 

Proposition 2.1. I/£ G T P TQ satisfies (12TTj) . thenTFL(^) G T FL{p) M l 
satisfies (12.31) . Therefore, if £ is a FLf-projectable solution of (12. ip . 
i/ien rfs projection TFLf(^) is a solution of (12. 3p . 

Conversely, if Y is a solution of (12. ip . t/ien an?/ -FL/ projectable 
vector field on Pf which projects on Y, is a solution of (12. 3p . 

Next, we shall discuss the SODE problem as it was stated by M.J. 
Gotay and J.M. Nester [HE]. 

The results can be summarized in the following result. 
Theorem 2.2. 

(i) If £ is a FLf-projectable vector field on Pf then for any p G M/ 
i/jere exists a unique point in each fiber FLJ 1 ^), denoted by 
Vt{p) a t which £ is a SODE. The point rj^{p) is given by 



is a section of FLf : Pf — > Mf and on Im(f3^) there exists a 
unique vector field, denoted by Y^, which simultaneously satis- 
fies the equations 



(ii) The map 




i Y .n L = dE L} SY^ = A. 
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We will now recall the construction of the solution of the dynamical 
equation which simultaneously satisfies the SODE condition. If Y — 
(F L f) , then Y is a vector field on Mf satisfying iy = dh\. The 
vector field Y^ described in (ii) is given by 

Y^(p)) = T^(Y(p)), for all p G M f . 
A detailed discussion can be found in [191 E El El fTTj . 



3. The Skinner and Rusk formalism 



Skinner and Rusk, [2^1 125] , have considered a geometrized framework 
where the velocities and the momenta are independent coordinates. 
Indeed, they considered the dynamics on the Withney sum of TQ (the 
space of velocities) and T*Q (the phase space). 

In this section we will briefly recall the Skinner and Rusk formalism. 

Let Q be a different iable manifold and L : TQ — > R a lagrangian. We 
can consider the bundle TQ @T*Q given by the Withney sum of tq : 
TQ ^Q and n Q : T*Q ->- Q. We will denote by pr\ ■ TQ®T*Q ->• TQ 
and pr 2 : TQ © T*Q — > T*Q the projections onto the first and second 
factors, and by pr : TQ © T*Q — > Q the projection onto Q. We then 
have the following commutative diagram 



TQ © T*Q 



(3.1) 




We can define a function 

T> : TQ © T*Q — ^ K 

a p ) ->> T)(X P , cKp) = a p (X p ) - L(X P 

In bundle coordinates (q A , v a ,pa), D is given by D(q A , v a ,pa) 
L(q A ,v A ). The function D is sometimes refered as the Pontryagin 
hamiltonian or generalized energy (see [26J). 

We can define a 2-form Q on TQ © T*Q by Q = pr^fig), where Qq 
denotes the canonical symplectic 2-form of T*Q. 

Next, we can consider the presymplectic system given by (Wo = 
TQ © T*Q, Q, dD) and study the equation 

i x tt = dD, (3.2) 

applying the Gotay-Nester-Hinds algorithm of constraints. Hence, we 
obtain 

W\ = {i£ W such that there exists X G T X W satisfying i x ^ = dD}. 



8 M. DE LEON, D. MARTIN DE DIEGO, AND M. VAQUERO 

In canonical coordinates (q A , v a ,pa), we have 
fi = dq A A dpA, 

■§^dq A ' v/'.i 777^ 



dD = -§kdq A + (pa ~ §k)dv A + 



So, given a tangent vector X = a A -^+b A -^+c A -^- e T^a^Wo 
we deduce that 

i x n = -c A dq A + a A dp A 
and (13 .2p is equivalent to the following conditions 

a A = v A , 

c A = (3-3) 
_ |^ = o, 1 < A < n. 

Next, we should restrict the dynamics to W\ = {{q ,v ,pa) 6 
W such that p A = J^}, that is, W\ = graph(FL), where FL :TQ ->■ 
T*Q has been defined in section |2j 

Accordingly with the Gotay-Nester-Hinds algorithm, a solution X 
must be tangent to W\. Assume that such X has the local expression 

Then, taking into account (13. 3p and (I3.4p .we deduce 
a A = v A 

r y ~r a„,APt„,B u 



dv A dq B dv A dv B dq A ' 

If there exists such a vector field X tangent to W±, satisfying the 
above conditions, we have done, and the final constraint manifold Wf 
is just W\. For instance, if the lagrangian is regular, det( d ®BQ v A ) ^ 0, 

—A 

we can compute b explicitly. If we denote by Cab the matrix Cab = 



(y d ®BQ v A ^j and C AB its inverse, then 



t = -C AB (v A ° 2L 01 



dv B dq A dq A J 

Otherwise, we need to continue the process, and then we obtain a 
sequence of submanifolds 

. . . ^ W k . . . W 2 Wi Wo = TQ © T*Q. 

If the algorithm stabilizes, that is, there exists k such that Wk = 
Wk+i, then Wu is called the final constraint submanifold and denoted 
by W f . 
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4. A Hamilton- Jacobi theory in the Skinner-Rusk setting 

In this section we will develop a Hamilton- Jacobi theory in the 
Skinner-Rusk formalism. We will use the same notation introduced 
in the previous sections and discuss separately the regular and the sin- 
gular cases. 

4.1. The regular case. Assume that we begin with a regular la- 
grangian L : TQ ->■ R. Then, W f = W v 

A section of TQ © T*Q is given by a = (Z, 7) where Z and 7 are a 
vector field and a 1-form on Q, respectively. Assume that a satisfies 
the following conditions 

(i) Im(cr) C W 1 = graph(FL), and 

(ii) d(pr 2 00) = dj = 0. 

Then, by the regularity of L, we know that there exists a unique 
vector field on W±, say X, satisfying 

%x ^ = dD, 

and then we can define a vector field on Q by 

X a (p) = Tpr{X(a(p))), for all p E Q. 

Now we have the following proposition. 

Proposition 4.1. Under the previous conditions, d(D o a) = if and 

only if the vector fields X and X a are o-related. 

Proof. 

(t s 11 

Assume that d(D o a) = holds, then we will prove first that 

(i(X-Ta(X°))Q = 0)|l m ( CT ). 

It is clear that if x e Im(cr) then T X {TQ@T*Q) = T x lm(a)+V, where 
V denotes the vertical bundle of the projection pr : TQ © T*Q — > Q. 
We will show that Z(x-t<t(x ct ))^ anihilates T x lm(o~) and V. Indeed, by 
the definition of Q, it is obvious that Q vanishes acting on two elements 
of V. Since X — Ta(X a ) is vertical, we have 

(i(X-Ta(X«))ty (V) = 0. 

Given p e Q, since X is a solution on Wi, we get 



(i x(p) n) o T<j(p) = TD(<j(p)) o T<j(p) = T(D o a)(p). 
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On the other hand, {iTa(x a ( P )) ° Ta{p) = since for any Y e T P Q 
we have 

(iTa(x^ p)) n) (Ta(p)(Y)) = n(Ta(X°(p)),Ta(Y)) 

= n(Ta(X°(p)),Ta(Y))=pr* 2 (n Q )(Ta(X°(p)),Ta(Y)) 

= (n Q )(Tpr 2 oTa(X°(p)),Tpr 2 oTa(Y)) = (Q Q )(T 7 (X-(p)), T 7 (y)) 

= -d 1 (T 1 (X°(p)),T 1 (Y)) 

= 

and so, we conclude that 

{i{x-T*(x°))ti) (TIm(a)) |Im(a) = 0, 

which implies 

(i {X -T* f (x°))ty (V + Tlm(a)) ]Im{a} 

= (i {X -Ta f{ x^) (T(TQ®T*Q)) llm{a) = 0. 

Therefore (X — Ta(X a )) e ker(f2). This means that i{x-Ta{x°)) ^ = 0, 
and hence t* (i(x-Ta(x°))Q) = i{x-Ta(x°)) = 0, where i : W\ — >■ 
Wo is the inclusion. 

It is not hard to see, that if L is regular then i*Q is symplectic and 
so(X = Ta(X-)) |Im(CT) . 

"<^" Since {(i {X -T*(x°))ty oTo = d(D oa)),HX = Ta(X CT ), then 
d(Doa) = 0. □ 

4.2. The singular case. Assume now that L : TQ — > R is an almost 
regular singular lagrangian. 

Suppose that the algorithm of Gotay-Nester-Hinds applied to (Wo = 
TQ(BT*Q, f2, dD) stabilizes at a final constraint submanifold W/. By 
construction, there exists at least one vector field X on Wf such that 

(i x ^ = dD)\ Wf 

We need some regularity conditions, thus we will also assume that 
Qi = pr(Wi) are submanifolds and that pri = pr\ w . : Wi — >■ are 
submersions. 

A section of pr : TQ © T*Q — ?• Q is given by a = (Z, 7), where Z 
and 7 are respectively a vector field and a 1-form on Q. We will denote 
by Of the restriction of a to Q/ = pr(W/) of er. Suppose that a verifies 
the following conditions: 

(i) Im(o-) C Wl 

(ii) lm(a f ) C W/. 

(hi) d(pr 2 o o") = (I7 = 0, that is, 7 is closed. 
Using a we can define a vector field on Qf by 

X"(p) = Tpr(X((7 / ( P ))), P6Q/. 
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The construction is illustrated in the following diagram 



W - W f — ^ TW f 

Tpr f 




The relation between T<jf(X a ) and X is shown in the following theo- 
rem. 

Proposition 4.2. The conditions 

d(Doa) lQf = 

and 

(X-Ta f (Xneker(Q)) lIm{(Tf) 

are equivalent. 

Proof. The proof follows by similar arguments as in Proposition 14. 11 

□ 

Definition 4.3. A section a ofTQ © T*Q, a = (Z,j), satisfying the 
following conditions 

(i) Im(a) C W x . 

(ii) Im(a f ) C W f . 

(iii) d{pr 2 o a) = dj = 0. 

(iv) d(Doa)\ Qf =0 

will be called a solution of the Hamilton- J acobi problem for the 
lagrangian L in the Skinner-Rusk setting. 

Remark 4.4. The last proposition says that Taf(X a ) is a vector field 
along Im(cry) which is also a solution of the equation (13. 2p . So if we 
find an integral curve c(t) of X a on Qf, then (a/ o c)(t) is an integral 
curve of a solution of (13.21) . o 

Remark 4.5. The natural question is if X and X a are <j/-related in the 
singular it happens in the standard Hamilton- Jacobi theory, 

see [T7j. The answer is that, as we discussed later (section [6]) , in some 
cases the fields are not necessarily <7/-related. o 

5. Comparison with the Hamiltonian and lagrangian 

settings 

In the previous section we have developed a Hamilton- Jacobi the- 
ory in the Skinner-Rusk setting. The Skinner-Rusk formalism unifies 
lagrangian and hamiltonian formalisms, so we would like to relate the 
present Hamilton- Jacobi theory to the corresponding ones for the two 
formalisms (see [I7J). 
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5.1. The hamiltonian setting. 

5.1.1. The regular case. If the lagrangian, L, is regular, that is, FL is 
a local diffeomorfism, then we can define locally a hamiltonian function 
h : T*Q — > R by h = E^oFL" 1 . Let us now assume that the lagrangian 
is hyperregular, that is, FL is a global diffeomorfism and h is globally 
defined. Denote by the corresponding hamiltonian vector field 

i Xh VL Q = dh. 

Let 7 be a closed 1-form on Q; then we can define a vector field on Q 
by 

X» = TttqIX^Ip))) for all p G Q. 

Then we have the following Hamilton- Jacobi theorem. 

Proposition 5.1. The vector fields X and X 1 are ^-related if and only 
if d(h o 7 ) = 0. 

Proof. For a proof see pQ . □ 

5.1.2. The singular case. Since we are considering an almost regular 
lagrangian L : TQ — > R, then we can apply the Dirac theory of con- 
straints developed in Section [2J 

We have to study the presymplectic system given by (Mi, Qi = 
j*^q, dhi), where j\ : Mi — > T*Q is the inclusion and hi is defined 
implicitly by hi o FL = E^. 

If we apply the Gotay-Nester-Hinds algorithm, we obtain a sequence 

. . . M k ^ ■ ■ ■ M 2 Mi T*Q; 

assume that we obtain a final constraint submanifold, denoted by My. 
We also assume that Qi = ■uq^Mi) are submanifolds and that 7Tj = 
7TQ| M : Mj — )■ are submersions. 

Remark 5.2. It is important to notice that the algorithm of Gotay- 
Nester-Hinds applied to the same lagrangian in the Skinner-Rusk set- 
ting and in the corresponding hamiltonian setting does not neces- 
sary stop at the same level. For example, the lagrangian given by 
L(q 1 , q 2 , v 1 , v 2 ) = v 1 q 2 produces the two presymplectic systems (Mi, Qi, 
dhi) and (W = TQ © T*Q,Q,dD). The first algorithm stabilizes in 
k — 1, but the second one does in k = 2. 

o 

Let 7 be a 1-form on Q satisfying the following conditions: 

(i) Im( 7 ) C Mi. 

(ii) Im(7/) C Mf, where 7/ denotes the restriction to Qf of 7. 

(iii) d~f = 0. 
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Then, if Y is a vector field on Mf solving the equation iy Qi = dh\, 
we can construct the vector field Y 7 on Q/ biven by 

Y\p) = T7T Q (Y( V (p))), for each p E Q f 

and obtain an analogous of theorem 14. 21 (notice that in this case we can 
ensure that the vector fields are 7/-related, see [TF] for the details). 

Proposition 5.3. We have 

d{h\ o j)\q. = Y and Y 1 are -related. 

Proof. Given q E Qf, we have 

(*(Y(7( g ))-T, 7/ (y*r( q ))) T <H = ^(7(9)) fl i T s7 - «T q7/ (y7( g )) fii o T 9 7 
= d lf{q) h 1 o T g 7 = o 7) 

where we have T g 7/(y 7 ) = ly^Y 7 ) and 

= ( 7 *fii)(^(g),F(g)) = rf 7 (^ 7 (g),^(g)) = 0, 

for all Y q E T 9 Q. 

The previous discussion can be applied to every point g G Q/; there- 
fore, taking into account that Vt\ vanishes acting on two vertical tangent 
vectors, we can deduce the following 

Y - r 7/ (YT) G ker(fi 1 ) ^ o 7 )| Q/ = 0. 

As we did before, we will see that Y and Y 7 are 7/ related. 
Remember that for any point p of M\ we have a decomposition 

T p (T*Q) = T p M 1 + V p (T*Q), 

where V(T*Q) denotes as above the space of vertical tangent vectors 
on p. 

Since Y — Tjf(Y rr ) is vertical at the points of Im( 7 j), given any 
U E V p , p E Im( 7 /), then 

n Q (Y-T 1 (Y^),U) = 

Now, given U E T P M\ we get 

Q Q (Y - r 7/ (Y"i% U) = Qt(Y - T 7/ (F 7 ), U) = 

because (F-T7 / (Y 7 )) E ker(fii), and hence tt Q (Y - T lf {Y~<), Z) = 
for any tangent vector Z G T P {T*Q) at any point of Im( 7 /). Since Oq 
is non-degenerate, we deduce that Y = T^fiY 1 ) along Im( 7 /). 

□ 

Definition 5.4. A 1-form 7 satisfying the previous conditions will be 
called a solution of the Hamilton- J acobi problem for L in the 
hamiltonian setting. 
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We are now going to relate the Hamilton- Jacobi problem in the 
Skinner-Rusk setting and the corresponding one in the hamiltonian 
setting. First, the following result gives the relation between W~i and 
Mi, and also a relation between solutions of equations ( 12. 3 j) and (13. 2|) . 

Lemma 5.5. 

(i) If X E T P W X satisfies i x ^ = dD, then X 2 = Tpr 2 {X) E 
T W2 (p)M\ satisfies i X2 & — dh\. 

(ii) For each step k of the constraint algorithms applied to the 
presymplectic systems (Mi, dh\) and (Wo = TQ © T*Q, 
Q, dD) we have 

pr 2 (W k ) C M k , 

and, if we denote the respective final constraint submanifolds 
by Wf and Mf, then 

pr 2 (W f ) = M f . 

(hi) We have pr{W f ) = n Q (M f ) = Q f . 

Proof. 

(i) Recall that a vector £ E T^a^a-jTQ, £ = u A -^ + w A -^ satisfies 
it il, = dE L iff 

9 2 L (v B_ u B ) = Q 



dv A dv B 

d 2 L )( _ 

dv A dv B dv A dq B dq A dv B dq A 



d 2 L B , d 2 L B _ dL_ _ d 2 L ( B _ ,.B\ 



If X E TpWi verifies ix ^ = dD, then X has the expression (13. 4p 
and satisfies (13.51) . So, it is clear that Xi = Tpr^X) satisfies i Xl &l = 
dE L . Since X is tangent to W t , X 2 = Tpr 2 (X) = TFL o Tpr x (X) = 
TFL(Xi) and using Proposition ^, ll we can conclude that ix 2 ^i = dh\. 

(ii) It will be proved by induction. 

For k = 1 we have pr 2 (W±) = Mi since W\ = graph(FL). 

Assume that pr 2 (W k ) C M k . Then 
Wk+i = {x E W k such that there exists X E T x Wk satisfying i x Q = dD} 
Mk+i = {y E Mk such that there exists Y E T y Mk satisfying iy £l\ = dh\}. 

If x E Wk+i, then there exists X E T x W k , satisfying i x Q = dD. 
Since pr 2 (W k ) C M k , Tpr 2 (X) E TM k and by (i) i Tp r 2 (X) ^i = dh x . 
Thus, we have proved that pr 2 (x) E M k+1 and that pr 2 (W k ) C M k . 

To prove that pr 2 (Wj) = Mf, take a solution Y of equation 12.31 on 
Mf. Then we can construct a vector field £ on Pf which is FLj-related 
with Y, and using Theorem 12.21 we obtain a vector field Y^ along the 
image of the section ^ which satisfies (12.11) and (12.21) . We can construct 
the map 

^: Mf — > TQ © T*Q 

(q A ,p A ) ->• (^(q A ,p A ),(Q A ,PA))- 
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It is easy to see, that the vector field T/3JY) on Im(/3^) is a solution 
of (13 .2p . By the maximality of the final constraint manifold Wf , we can 
conclude that Im(/^)c Wf, but M f = pr 2 (Im(^)) C pr 2 (W } ) C M f 
and then the result follows. 

(iii) It is a direct consequence of (ii) and the commutativity of dia- 
gram (13.11) . □ 

A solution of the Hamilton- Jacobi problem as stated in the previous 
section is given by a section a of TQ © T*Q, so a = (Z,j), where Z 
and 7 are a vector field and a 1-form on Q, respectively. 

We will see that 7 satisfies the Hamilton- Jacobi problem in the hamil- 
tonian sense. 

From the fact that a is a solution of the Hamilton- Jacobi problem 
in the Skinner- Rusk setting, we deduce: 

(i) Since Im(a) C Wi, then Im(7) = pr 2 (Im(a)) C pr 2 (Wi) = 

(ii) Since Im(<7/) C Wf, then Im(7/) = pr 2 (Im((Tf)) C pr2(Wf) 
= M f . 

(iii) Since d{pr 2 o a) = d^y = 0, then 7 is closed. 

(iv) Since Im(er) C W\, then D o a = hi o 7 and then, using that 
d(D o a)\ Qf = 0, we finally get d(h ± o j)^ Qf = 0. 

On the other hand, given a vector field X on Wf which is a solution 
of (13. 2p . we can obtain a solution of (I2.3P along Im(7/) by defining 

X 2 ( 7 /(p)) = Tpr 2 (X(a f (p))), for all p E Q f . 

Now, from Lemma 15.51 it follows that X 2 is a solution of (12.31) . 

As above we can construct the projected vector field on Qf, by 
putting 

XUp) = Tn f (X 2 ( V (p))), for all p G Q f . 

Remark 5.6. By the commutativity of the diagram (13. ip we deduce 
that pr = ttq o pr 2 , and in consequence we have 

X°(p) = Tpr(X(a f (p))) = Tn Q o pr 2 (X(a f (p))) = Tn f (X 2 ( V (p))) 

for all p 6 Qf, and so, X° = X 2 - o 

Summarizing the above discussion, we can conclude that it is possible 
to relate the Hamilton- Jacobi theory in the Skinner- Rusk setting to the 
Hamilton- Jacobi theory on T*Q. In this case the vector fields X 2 and 
XI are 7/-related. 

5.2. The lagrangian setting. In this section we will relate the Hamil- 
ton- Jacobi theory developed in the Skinner-Rusk setting with the cor- 
responding one on the lagrangian side 
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5.2.1. The regular case. If the lagrangian L is regular, then we have a 
symplectic system given by (TQ, Q L , E L ). Then there exists a unique 
solution £ of the equation 12.11 which automatically satisfies the SODE 
condition. 

Given Z a vector field on Q such that Z*Q^ = we can define the 
following vector field on Q 

e{p) = TT Q (Z(Z(p))) for all p E Q 

and obtain the following result. 

Proposition 5.7. Under the previous conditions, the vector fields £ 
and £ z are Z -related if and only if d(Ei o Z) = 0. 

Proof. The proof is a consequence of Proposition 15.11 □ 

5.2.2. The singular case. In this case, we will discuss the presymplec- 
tic system given by (TQ,£Il,(IEl). Applying the Gotay-Nester-Hinds 
algorithm we obtain a sequence of submanifolds 

• • • P k ^ ■ ■ ■ P 2 Pi = TQ. 

We also assume that Qi = tq^Pi) are submanifolds and that T{ = 
prQ| P : Pi — > Qi are submersions, for any index i. 

Remember that the algorithm of Gotay-Nester-Hinds applied to the 
presymplectic systems (M 1; Qi, dh\) and (TQ, Q L , dE L ) stop at the 
same step, so we will denote the final constraint manifold of the system 
{TQ,n L ,dE L ) by P f . 

Let Z be a vector field on Q satisfying the following properties: 

(i) lm(Zf) C Pf, where Zf denotes the restriction of Z to Qf. 

(ii) Z*Q L = 0. 

Then, if £ is a vector field on Pf solving the equation z^Ql — dEL, 
we can construct the vector field £ z on Qf by 

i z {p) = Tr Q (aZ f (p))), for all p E Qf. 

Now, we can develop the corresponding Hamilton- Jacobi theory in the 
lagrangian setting. 

Proposition 5.8. Under the above hypothesis for Z we have 
d(E L o Z)\ Qf = & (£ - TZf(i z )) E ker(fi L ). 

Proof. "=>" 

Assume that d(Ei o Z)\q. = holds, then we will prove that 

For any x E lm(Zf) we have the decomposition T X (TQ) = T x Im(Z) + 
V{TQ), where V(TQ) denotes the vertical bundle of the projection 
t q '■ TQ -> Q. 
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Since Q L vanishes acting on two elements of V(TQ) and £ — TZj 
is vertical, we have 

Since £ is a solution along Im(Zj), we have 

(i €(p ) Q £ ) o 7V(p) = d zip) E L o TZ/(jo) = <^,(.E L o Z) 
for any p E Qf. 

On the other hand, (i TZf (^z^ Q L ) o TZ{p) = 0, since for any F G 
T p Q we get 

(i( € _ rz , (€ *)) fi L ) oTZ(p)(Y) = n L (TZ f (£ z (p)),TZ(Y)) 
= (Z*n L )((Z z ,Y) = -*y(Z z ,Y) = 
and so we can conclude that 

{\i-Tz s (^)) n i) ( Tlm ( z )) = 

Since i( € _ rZ/K z))^L = o Z), if (f - TZ f (£ z )) e ker(fi L ), then 
d(E L oZ) lQf = □ 

Definition 5.9. yl vector field on Q, Z satsifying the previous condi- 
tions will be called a solution of the Hamilton- Jacobi problem for 
L in the lagrangian setting. 

The vector fields £ and £ z are not necessarily related as the next 
example shows. 

Example 5.10. Let L : TM 2 — > R be the lagrangian given by 
L(q l ,q 2 , v 1 , v 2 ) = q l v 2 + q 2 v l 

We have 

FL(q\q 2 ,v\v 2 ) = (q 1 , q 2 , q 2 , q 1 ), 

E L (q 1 , q 2 , v l ,v 2 ) = q 1 v 2 + q 2 v 1 — q 1 v 2 — q 2 v 1 = 0, 

^l = 0, 

so every vector field £ on TM 2 satifies 

i ( n L = dE L . 

Therefore, the algorithm of Gotay-Nester-Hinds stabilizes at the first 
step, and P f = P 1 = TQ. 

Moreover, every vector field Z on R 2 is a solution of the Hamilton- 
Jacobi problem, since E L o Z = and Z*Vl L = 0. 
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Let £ be the solution satisfying the SODE condition given by 

d d 

£(V, q 2 , v , v 2 ) = v 1 h v 2 1 1 

dq 1 dq 2 dv 1 dv 2 

Let Z be 

Z ^^ ) = W + dq- 2 



An easy computation shows that 



i 2 \ 9 d 



but 

Thus, the vector fields £ and £ z are not Z-related. 

Next we will show that a solution of the Hamilton- Jacobi problem in 
the Skinner-Rusk formalism induces a solution of the Hamilton- Jacobi 
theory in the lagrangian setting. 



The following lemma is analogous to Lemma 15.51 
Lemma 5.11. 

(i) If X e TpW\ satisfies %x ^ = dD, then X 2 = Tpr\{X) satisfies 
ix 1 &l = dEi and the SODE condition (12. 2p . 

(ii) For each step k of the constraint algorithm applied to the presym- 
plectic systems {M u Q. x , dhi) and (W = TQ®T*Q, Q, dD), 
we have 

pr 2 {W k ) C P k 

(iii) We have pr(Wf) = tq(P/) 

Proof. 

(i) and (ii) are proved using similar arguments to that in Lemma [531 
(iii) Since the following diagram 




is commutative, and FLf is a surjective submersion, we deduce that 
7r Q (M f ) = r Q (P f ). By Lemma |53J (iii), we obtain ir Q (M f ) = pr(W f ), 
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and the result follows. The situation can be summarized in the follow- 
ing commutative diagram 




□ 

If a = (Z, 7) is a solution of the Hamilton- Jacobi problem, we deduce 
the following results: 

(i) Since Ihi(<7/) C Wf, then pri(Im(cx/)) C pr\(Wf) C Pf. 

(ii) We have Z*Vt L = 0, since Z*Q L = Z*(d6 L ) = d{Z*9 L ) = 
d{FL{Z)) = d 1 = 0. 

(hi) Since Im((j) C Wi, then D o a{p) = El o Z(p) and, because 
d{D o a)\ Qs = 0, then d[E L o Z)\ Qf = 0. 

Now, given a solution X of (13. 2p . we can obtain a solution of (12. ip 
along lm(Zf) using Lemma [5.111 and putting 

X 1 (Z f (p)) = T pri (X(a(p))), for all p G Q f 

We can also define the vector field on Qf given by 

X?(p) = TT Q (X 1 (Z f (p))). 

The vector fields X\ and Xf are not ^/-related in general, as we 
have proved in example 15.101 

Remark 5.12. By the commutativity of diagram (13. ip we have pr 
= tq o pri and hence 

X°(p) = Tpr(X(a f (p))) = Ttq opn(X(Zf (p) ) ) = Tr Q {X x {Z f {p))) , 

for all p G Qf, and so X a = Xf . 

Moreover, since X\ satisfies the SODE condition, then 

*f(p) = TrqiX^Zip))) = ^(^(^(p))) = Z(p) = ^(p), 

and we have 

X CT = X% = XI = Z f . 

Note that this means that we only need to compute X^ to obtain 
Z f . o 
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6. Final considerations 

In the last section we show that a solution of the Hamilton- Jacob i 
problem in the Skinner- Rusk setting, a = (Z, 7), gives a solution of the 
Hamilton- Jacobi problem in the lagrangian and hamiltonian settings 
(Z and 7 respectively). A solution of the equation 13 . 2 1 along Im(cr) can 
be also projected to solutions of 12.11 and 12.31 along Im(Z) and Im(7), 
denoted respectively by X\ and X 2 . 

If we take a vector field X solution of the equation (13 .2p on Wf, using 
a we can compute X a . Now we can easily conclude that the vector 
fields X and X 1 are 07 related iff the corresponding vector fields Xi 
and Xf are Zf related in the lagrangian setting. 

To illustrate the above results we revisite example 15. 101 in the Skinner- 
Rusk setting and apply the corresponding Hamilton- Jacobi theory. 

Example 6.1. Consider the lagrangian given in Example 15.101 

L(g\ q 2 , v 1 , v 2 ) = q 1 v 2 + q 2 V 1 . 

Then, on TW 2 © T*R 2 we have 

D(q 1 , q 2 ,v 1 ,v 2 ,p l} p 2 ) = v 1 p 1 + v 2 p 2 + v l q 2 + v 2 q 1 , 

and hence 

dD^q 1 , q 2 , v 1 ,v 2 ,pi,p 2 ) = —v 2 dq 1 — v 1 dq 2 + {jp\ — q 2 )dv Y 

+(Pi — q^dv 2 + v 1 dpi + v 2 dp 2 

Recall that we must compute 
Wi = {(g A , v A ,p A ) such that there exists X G T (9 a^a w) TI 2 © T*R 2 
satisfying ix^l = dD}. 

If 

^ a dq 1 a dq 2 ^ dv 1 ^ dv 2 ° dpi ° dp 2 

then 

ix ^ = —c^dq 1 — c 2 dq 2 + a l dp\ + a 2 dp 2 (6.3) 

and so 

a 1 = v 1 , a 2 = v 2 , c 1 = v 2 , c 2 = v 1 , p\ — q 2 = 0, p 2 — q 1 = 

(6.4) 

must hold. 

Therefore, W\ = {(q l , q 2 , v l , v 2 , q 2 , q l ) such that q A , v A eR} = graph(FL). 
Next, we compute 

W 2 = {(q 1 , q 2 , v 1 , v 2 , q 2 , q 1 ) G W\ such that there exists 
X G TLi^y ^a^gi^Wi satisfying i x ^ = dD} 
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If X e TWi then X can be locally expressed as 

I / d 2 L 1 , d 2 L Q 2 9 2 L , d 2 L \ a 
^cfe 1 ^ 1 dv x dq 2 dv 1 dv 1 dv 2 dv 1 ' dpi 

^dv 2 dq 1 dv 2 dq 2 dv 2 dv 1 dv 2 dv 2 ' dp2 

- 1 JL + i hl_d_ + A2^_ , fl 2 j9_ _j_ 1_9_ 



(6.5) 



Taking into account (16. 4ft and (16.51) . for every point (q 1 , q 2 , v 1 , v 2 , g 2 , g 1 ) 
G Wi we obtain 

X - v 1 ^ + w 2 ^ + 6 1 ^ + 6 2 ^ + v 2 ® + u 1 ^ 
5g 2 dg 2 df 1 dv 2 (9pi (9^2 

for arbitrary 6 1 , b 2 , and so = W\ and therefore the final constraint 
submanifold is W\\ consequently, Qf = Q. 

Now, a solution of the Hamilton- Jacobi problem in the Skinner-Rusk 
setting is given by a = (Z, 7) such that 

(i) Im(cr) C Wx. 

(ii) Im((j j) C Wf. 

(iii) d{pr 2 o a) = = 0, that is, 7 is closed. 

(iv) d(Doa) ]Qf = 

It is easy to see that every pair given by a vector field Z and its 
image by the Legendre transformation, that is (Z, 7 = FL(Z)) is a 
solution of the problem. In fact, by construction Im(a") C W\ and 
D\Wi = 0^Doa = 0. Following the argument in example 15.101 we 
can take Z(q 1 , q 2 ) = + and so 

d d 

If we consider the solution 

d d d d d d 

X{q\ q 2 , v\ v 2 ) =v 1 — + v 2 — + — + — + v 2 — + V 1 — 
oq L oq z ov 1 ov z opi op2 



then 



and 



Ta{X q ,q ^ + ^ + ^ + 



A direct inspection shows that 
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We can also obtain information of the Hamilton- Jacobi problem in 
the Skinner-Rusk setting from a solution of the Hamilton- Jacobi prob- 
lem in the hamiltonian side. 

If 7 is a solution of the Hamilton- Jacobi problem in the hamiltonian 
setting and Y a vector field on Mf wich is a solution of equation (12. 3p , 
then we can define Y 1 as before. 

We can also define a section a of prf : Wf —> Qf given by a(p) = 
(Y" / (p), r y{p)) for all p G Qf. An easy computation shows that Tcr(F 7 ) 
is a vector field along Im(a) which solves (13. 2p . Moreover if we find a 
vector field Z on Q such that FL o Z = 7 and Zf = Y 1 , then the pair 
(Z, 7) is a solution of the Hamilton- Jacobi problem in the Skinner- Rusk 
setting. 

7. Appendix: The Gotay-Nester- Hinds algorithm of 

constraints 

In this section we will briefly review the constraint algorithm of con- 
straints for presymplectic systems (see [TT1 17]). 

Let Mi be a manifold, Q a presymplectic structure on Mi, i.e., Q 
is a closed 2- form, and a a 1-form on Mi. We will call (Mi, £1, a) a 
presymplectic system. 

Gotay et al. developed an algorithm to find N, a submanifold of Mi 
where we can solve the equation 

i x tt = a (7.1) 

with X tangent to N. 

The previous equation could not hold for every point of Mi , because 
a could not be in the range of Q. So it is necesary to introduce the 
following set 

M 2 = {p E Mi such that there exists X e T p Mi satisfying ix f2 = a}, 
and it is assumed that M2 is a submanifold. 

At the points of M 2 there exists solution to equation (17. II) but in an 
algebraic sense, that is, the solution could not be tangent to M 2 . This 
forces a further restriction to 

M3 = {pG M 2 such that there exists X G T p M 2 satisfying i x f2 = a}, 
which is also assumed to be a submanifold. 

Proceeding as above, the algorithm will produce a sequence of sub- 
manifolds 

• • • M 3 . . . -V 3 M 2 -V 2 Mi 

where 

M/ + i = {p G M[ such that there exists X G T p M t satisfying i x £1 = a}, 
and ji denote the inclusions. 
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There are three possibilities: 

(i) There exists k such that = 0. 

(ii) There exists k such that = Mk+\. 

(iii) The algorithm does not end. 

In the second case the submanifold M/, is called the final constraint 
submanifold and is denoted by Mf. By construction there exists a 
vector field on Mf such that is solution of equation (17.11) . The third 
case is only possible in the infinite dimensional setting. In this case, 
the final constraint submanifold is defined by Mf = flj = iMj. 

Note that the final constraint submanifold is maximal in the sense 
that if R is submanifold of Mi where there exists a tangent solution of 
equation (17. lj) . then R C Mf. 
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